Integral equation for inhomogeneous condensed bosons generalizing the 
Gross-Pitaevskii differential equation 
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We give here the derivation of a Gross-Pitaevskii-type equation for inhomogeneous condensed 
bosons. Instead of the original Gross-Pitaevskii differential equation, we obtain an integral equation 
that implies less restrictive assumptions than are made in the very recent study of Pieri and Strinati 
[Phys. Rev. Lett. 91, 030401 (2003)]. In particular, the Thomas-Fermi approximation and the 
restriction to small spatial variations of the order parameter invoked in their study are avoided. 



PACS numbers: 03.75.Ss, 03.75.Hh, 05.30.Jp 

In a very stimulating recent contribution, Pieri and 
Strinati (referred to as PS below) P] have 'derived' 
the non-linear Gross-Pitaevskii differential equation for 
condensed bosons by taking as their starting point the 
Bogoliubov-de Gennes equation for superfluid fermions. 

The purpose of this Brief Report is to demonstrate 
that one can generalize the zero-temperature differen- 
tial Gross-Pitaevskii equation while remaining within the 
original framework of PS, an integral equation formula- 
tion then resulting. The framework of PS is provided by 
the coupled integral equations involving Green functions 
G21, Gil and Go. The equations are: 

Gii(r,r';L^,) = Go(r,r';^,)+ j dr"Go(r, r"; w,) 

xA(r")G2i(r",r';c^,), (la) 

G2i(r,r';w,) = - j dv" Go{v" ,v- -lo^) 

xA*(r")Gii(r",r';c.,), (lb) 

where u)s = (2s + l)7r//3 (s is an integer) is a fermionic 
Matsubara frequency, (3 — 1/k^T, Gn is the normal and 
G21 is the anomalous single-particle Green function. The 
third Green function appearing in Eqs. namely Go, 
satisfies the equation 

[luj, - H{v)]Go{v, r'; lo^) - <5(r - r'), (2) 

where the single-particle Hamiltonian H{v) is defined by: 



(3) 



/i being the Fermionic chemical potential. As PS stress, 
Eqs. (P), when taken together with the self-consistency 
equation for the gap function: 



are entirely equivalent to the Bogoliubov-de Gennes 
equations that describe the behavior of superfluid 
fermions in the presence of an external potential. Equa- 
tions (^EJ define what we have termed the original 
framework of the PS study. The constant Vq < enter- 
ing Eq. Q arises from the contact potential Vq5{y — r') 
assumed by PS to act between fermions with opposite 
spins. We also retain here their use of the ratio A(r)//i 
as an expansion parameter which allows the rapid trun- 
cation of such series, which then leads for strong coupling 
to an integral equation for the gap function 



^A*(r) = 



j drig(r,ri)A*(ri) 

-I- y dridr2dr3 i?(r,ri,r2,r3) 
xA*(ri)A(r2)A*(r3), 



(5) 



A*(r) 



G2i(r,r;ws 



(4) 



where R is written explicitly in terms of Go(r, ri;^^) 
in Eq. (15) of PS. However, as will emerge below, it is 
the non-local kernel Q{v, r') which is at the heart of the 
present study. In terms of the Green function Go entering 
Eq. ©, g(r,r') is given by [PS: Eq. (14)]: 

Q(r, r') = ^ E ^°(^'' ^' -^s)Go{v', r; u,). (6) 

We take the integral equation ISJ for the gap function as 
the starting point of this Brief Report. For our purposes 
below, it is then crucial to gain insight into the kernel Q 
in Eq. |(HJ , and in particular to carry out the summation 
explicitly over the Matsubara frequencies los ■ 

To gain orientation, let us first perform this summation 
when the external potential T^(r) is set to zero in Eq. ||2J). 
Having achieved this summation, we shall present a gen- 
eral method to allow the sum over Ug to be achieved for 
I^(r) 7^ 0, using earlier work of Stoddart, Hilton and 
March [3. 
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FIG. 1: Solid lines show r'^Qo{r'), where Q{r) is defined by 
Eq. @, as a function of r' — kpr, for several temperatures, 
given by /?' — Pjj, = 10 — 30 (bottom to top). Dashed lines 
are the asymptotic expansion Eq. Units are such that 

fc|/(4/i) = 1. 



Returning to the explicit form of Q{r,ri) given in 
Eq. JBJ above, it is natural to study first the transla- 
tional invariant, free-electron limit of Eq. ©, say Qo{i'), 
with r = |r — ri|, which is obtained by 'switching off' the 
one-body potential V{y). This amounts to replacing Go 
in Eq. © with the free-electron Green function Go. For 
the Fourier transform of Qoir), we formally find 



go(fc) 



dk' 1 - nF(^k-k') - "■F(fk' 



(27r)3 
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k-k' 



^k 



(7) 



where = k'^ /2m — jj, and nF(C) is the Fermi-Dirac dis- 
tribution function. However, it should be noted that, in 
three dimensions, Eq. contains a divergent contribu- 
tion at large wave-numbers, which implies a divergent 
behavior of Qo{r) at small distances r. Indeed, we find 
the asymptotic expansion (see also Ref. Q): 



1 1 



4^2 j.l2pi gij^hfj' 



» 1, 



(8) 



where r' = k-pr, kp is the Fermi wave-number, defined 
by ^ = fc|/2m, /?' = /3/i, and a = r'n/P'. Fig. ^shows 
then our numerical results for r"^Qo{r'), as a function of 
r', for several temperatures (/3' = 10 — 30). 

Following Stoddart et al. 0, the canonical density 
matrix G(r, r', /3) is defined by 



G(r,r',/3)= VV,(r)V'*(r')e-^^' 



(9) 



where /3 = l/k^T. Within the perturbative approach of 
March and Murray 0,01 with plane waves as the unper- 
turbed solution, the canonical density matrix can then 
be written to all orders in the external potential V{r) in 
terms of the free-particle canonical density matrix given 

by 



Go(0,/?) = (27r/3)-3/2exp(-zV2/3), 



(10) 



C(r,ro,/3) 



dzzCo{z,P)f{z,r,ro), 



(11) 



where / satisfies the integral equation 01 : 



fiz,r,ro) 



1 



S{z 



rol) 



^(ri) 



z ' ' J "27r|r — ri| 

x/(z-|r-ri|,ri,ro). (12) 



The desired Green function Go is then to be obtained 
from / entering Eqs. ((TT|l and itT^ as 



Go(r,ri;fc) 



dz z Go{z; k)f{z, r, ri 



where 



Go{z;k) 



Attz 



(13) 



(14) 



One may also take advantage of the expression in Eq. H13() 
of Go in terms of Go to rewrite the kernel Q{r, ri) defined 
by Eq. © as 

(5(r,ri)= / dzidz2ZiZ2f{zi,ri,r)f{z2,ri,r)Qo{zi,Z2), 
Jo 

(15) 



where the Fourier transform of Qo(zi, Z2) is given by 
1 - nF(^ki) - "-pCCks) 



(5o(ki,k2) = 



Cki + Ck2 



(16) 



Hence, the sum over Matsubara frequencies has still been 
carried out in the presence of an external potential ^(r) 
entering Eq. (|12|l for the function /. 

Because of current interest in harmonic confinement in 
magnetic traps at low temperatures, let us illustrate the 
rather formal Eqs. 111|) and (|12|l when the external po- 
tential V{r) has the explicit isotropic harmonic oscillator 
form in three dimensions, namely 



y(r) 



1 2 2 

—muj r . 
2 



(17) 



Following the pioneering work of Sondheimer and Wil- 
son (6| on free electrons in a magnetic field, the diagonal 
element G(r,r,/3) when V{r) is given by Eq. (|17(l takes 
the form (see e.g. 0, P- 27; see also 0) 



3/2 



sinh huj/S 



3/2 



X exp 



TO 2 ■ 1 



^^ujr^ tunh -hujl3 ] , (18) 



which is the so-called Slater sum of quantum chemistry 
(Fig.©. 

From Eqs. (|10|) and Ijlll) . performing the substitution 
t — z^/2, it then follows that f{z, r, rp) can be expressed 
as the inverse Laplace transform 



/(z,r,ro) = (2^)3/2/:-i \s-^/^C{r,ro, s'') 



(19) 
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FIG. 2: Shows diagonal element of the canonical density ma- 
trix C(r, r,/9), Eq. 11811 . and its Thomas- Fermi approxima- 
tion, Eq. H2()|l . as a function of r, for /3 = 0.5. Energies are in 
units of hiu, while lengths are in units of {h/mui)^^^ . 
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FIG. 3; Shows diagonal f{z, r, r) corresponding to the har- 
monic potential, as given by the inverse Laplace transform, 
Eq. 11911 . as well as the regular part of its Thomas- Fermi ap- 
proximation, Eq. 12111 . as a function of r, for fixed z — 0.5. 
Units as in Fig. ^ 



where (t, s) are conjugate variables with respect to the 
Laplace transform. 

Within the Thomas-Fcrmi (TF) approximation, we 
take: 



CTF(r,r,/3) = 



1 



(27r/3)3/2 



exp[-(3V{r)], (20) 



which is plotted also in Fig.|21for V^(r) given by Eq. (|17|l . 
For the value of /3 shown, the TF form Eq. H2U|) is seen to 
be a useful approximation to the exact result, Eq. (|18|l . 
Inserting Eq. (^0)1 into Eq. ifT^ we find 



/TF(.,r,r) = ^ - ^^MVm^z], (21) 

z z 

where Ji denotes the Bessel function of the first kind and 
order one. 

Fig- Elshows f{z, r, r) as a function of r for fixed z, as 
obtained by numerically performing the inverse Laplace 



transform in Eq. (|19|l for the harmonic potential case. 
The regular contribution to the analytic result for the 
Thomas- Fermi approximation, Eq. H21|l. is also plotted 
for comparison. The similarity in shape between approx- 
imate and (numerically) exact results for this harmonic 
confinement model seems to us rather remarkable. After 
this model test of a TF-like approximation invoked by 
PS, we return to the general case, based on the exact 
result Eq. (|15|l for the kernel Q{r,ri). 

Then, we invert the argument of PS but still use a 
further essential assumption of their study, namely that 
the condensate wave function $(r) entering the Gross- 
Pitaevskii equation is related to the gap function A(r) 

by 



$(r) 



Stt 



1/2 



A(r) = A:A(r). 



(22) 



Here, in the strong coupling limit, and following PS, ap ~ 
(2m|/i|)~^/^ represents the characteristic length scale for 
the non- interacting Green function Go , equal to Go above 
when V{r) is put equal to zero. 

Given the validity of this PS assumption, Eq. I22|l . we 
then rewrite Eq. © as an equation for $(r): 



Vo 



$*(r) 



J drig(r,ri)$*(ri) 

^ J dYidY2dYzRir, ri,r2,r3) 



fc2 

X$*(ri)$(r2)$*(r3). 



(23) 



This then is the proposed generalization of the Gross- 
Pitaevskii equation, but with (5(r,ri) to be calculated 
more accurately than by the Thomas-Fermi-like assump- 
tion of Fieri and Strinati l], via Eqs. Ijl5l) and H12|l . 

While Eq. (|23|l is a direct consequence of the above 
arguments, it remains an expansion in <i>, in suitable re- 
duced form. Therefore, a first attempt to simplify this 
Eq. H23|) is to retain the approximation given by the Pieri- 
Strinati approach in the 'smallest' term involving 0($'^) 
on the right-hand side of the basic Eq. (|23l) . Thus one 
reaches the (still non-local) equation for the condensate 
wave function <i>(r): 



--l<f*(r) = j driQ(r,ri)$(ri) 



-^|$(r)pa>(r) 



(24) 



For sufficiently small spatial variations in the conden- 
sate wave function <i>(r) in Eq. (|24|l . the basic nonlocality 
can be removed by Taylor expanding $(ri) around the 
position r in the integral term. This then characterizes 
the problem in terms of 'partial moments' of the kernel 
(3(r,ri), namely /(5(r,ri)dri and / (5(r, ri)|r - ri pdri. 
Such partial moments then enter the original Gross- 
Pitaevskii equation, as stressed by PS. 

In summary, we propose the retention of the non-local 
kernel Q(r,ri) as in Eq. (|23|l above, since the sum over 
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Matsubara frequencies in Eq. (jSj) has been performed in 
Eq. (|15l) . which is a central result of the present study. 
However, in the terms of 0{^^) in Eq. (|23|l . a sensible 
starting point is to follow the PS approximation displayed 
in Eq. 

As to future directions, evaluation of the non-local ker- 
nel in Eq. H15|l for other external potentials than the har- 
monic case in Eq. (|17|l is of obvious interest. For this 
latter model, though our Fig. 13 considers the diagonal 
element off{z, r, ri), the off-diagonal form of C(r, ri, /3) 
is known and numerical Laplace inversion to obtain 
/(z,r,ri) is entirely feasible. Then Q{r,ri) can be ob- 
tained, though of course numerically. 

The Gross-Pitaevskii equation is valid in the strong- 
coupling limit of superfluidity. It has to be stressed 
that in the weak coupling limit one can also de- 
rive a Ginzburg-Landau equation starting from the 
Bogoliubov-de Gennes equations. We note specifically in 
this context that the derivation of the Ginzburg-Landau 
equation in the weak-coupling limit for the harmonic trap 
was presented by Baranov and Petrov The results 



presented in this Brief Report are also relevant to the 
weak-coupling limit of superfluidity. 

Finally, we should mention the very recent discussions 
of the foundations of the Gross-Pitaevskii equation by 
Leggett ilflj. He concludes that there is no correlated 
many-body wave- function underlying their original equa- 
tion. It will be interesting for the future to know whether 
the non-local versions of Eqs. H23|l and (|24(l proposed here 
are still subject to this limitation. 
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